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A numerical study of the transfer across random fractal 
surfaces shows that their response are very close to the re- 
sponse of deterministic model geometries with the same frac- 
tal dimension. The simulations of several interfaces with pre- 
fractal geometries show that, within very good approxima- 
tion, the flux depends only on a few characteristic features of 
the interface geometry: the lower and higher cut-offs and the 
fractal dimension. Although the active zones are different for 
different geometries, the electrode responses are very nearly 
the same. In that sense, the fractal dimension is the essential 
"universal" exponent which determines the net transfer. 
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Many random processes such as aggregation, diffusion, 
fracture and percolation, build fractal objects Frac- 
tal geometry essentially describes hierarchical structures 
[H . If properties of these random systems depend on the 
hierarchical character of their geometry, then the study of 
a deterministic structure with the same fractal dimension 
may provide a good approximation of the random system 
properties [Q. The question is significant since fractal 
and pre-fractal geometries are widely used in mathemat- 
ical approaches or numerical simulations as a convenient 
model of irregularity. They are also more simply ad- 
dressed by algebraic calculations and incorporated into 
numerical models for computer simulation. It is then an 
important matter to decide whether simple determinis- 
tic, artificial, fractals could help determine the proper- 
ties of random, natural, fractals |^^. In particular, it 
is a question whether experiments performed on model 
fractal geometries may help understand the behavior 
of real complex structures. 

The property which is discussed here is the Laplacian 
transport to and across irregular and fractal interfaces. 
Such transport phenomena are often encountered in na- 
ture or in technical processes: properties of rough elec- 
trodes in electrochemistry, steady-state diffusion towards 
irregular membranes in physiological processes, the Eley- 
Rideal mechanism in heterogeneous catalysis in porous 
catalysts, and in NMR relaxation in porous media. In 
each of these examples, the interface presents a finite 
transfer rate, like a redox reaction, or a finite permeabil- 
ity, or reaction rate which is due to specific physical or 
chemical processes. 

The mathematical formulation of the problem is sim- 
ple. One considers the current flowing through an elec- 



trochemical cell as shown in Fig. ^ The current J is 
proportional to the Laplacian field VV, which can be 
viewed as an electrostatic field in electrochemistry, or a 
particle concentration field in diffusion problems. Then 
the flux and field are related by classical equations of the 
type J = — (jVF, where a is the electrolyte conductiv- 
ity (or particle diffusivity in diffusion or heterogeneous 
catalysis). The conservation of this current throughout 
the bulk yields the Laplace equation for the potential V\ 
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FIG. 1. Schematic representation of an electrochemical cell. 

The boundary presents a finite resistance to the cur- 
rent flow. In the simplest case, this resistance can be 
expressed by a linear relation linking the current density 
across the boundary to the potential drop across that 
boundary. The local flux and potential drop are then 
linked by transport coefficients, like the faradaic resis- 
tance in electrochemistry, the membrane permeability in 
physiological processes, or again the surface reactivity in 
catalysis. If one assumes that the outside of the irregular 
boundary is at zero potential, current conservation at the 
boundary leads to the following relation: 
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The parameter A is homogeneous to a length. Given 
the geometry, the value of this parameter determines the 
behavior of the system [||,^. The overall response of 
such a system is measured by one scalar quantity, its 
impedance Ztot^ which is the ratio between the applied 
potential and the total flux : 
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The contribution of the finite interface resistivity to 
this global impedance is given by a "spectroscopic" 
impedance, defined as: Z 



spec. 



Ztot — Zq, Zq be- 
ing the impedance of the cell with zero interface resis- 
tivity The main result discussed below is that the 
electrode impedance Zspec.^ is nearly independent of the 
random character of the fractal interface, even though 
the regions where the current is concentrated are very 
different. This is found from a numerical comparison 
between impedances of deterministic and random elec- 
trodes with the same fractal dimension. Two cases are 
studied: (a) deterministic and random von Koch elec- 
trodes (dimension Df = In4/ln3), (b) a deterministic 
electrode of dimension Df = 4/3 and a self- avoiding ran- 
dom walk geometry with the same dimension. 

The deterministic von Koch curve, or classical 
snowflake curve, is obtained by dividing a line segment 
in three equal parts, removing the central segment and 
replacing it by two other identical segments which form 
an equilateral triangle. A random von Koch curve can 
be defined simply by choosing randomly the side of the 
segment where the triangle is created at each step of the 
building process. This is shown in Fig. ^. After three 
or more generations, it looks more like a realistic ran- 
dom boundary than a simple mathematical curve. It is 
then possible to automatically generate different bound- 
aries that have the same fractal dimension and the same 
perimeter. By definition fractal geometries exhibit a 
large scale of lengths. For instance, at the sixth gen- 
eration, the ratio between the smallest feature / (smaller 
cut-off) of the irregular boundary and the diameter L 
(larger cut-off) is 3^ = 729 while the length of the 
perimeter is Lp = 4^1 = 4096/. Computing on a regu- 
lar grid within such geometries would be very memory 
and time-consuming. A finite element method is then 
used. The standard variational formulation of the prob- 
lem is discretized with a triangular mesh, obtained from 
a Delaunay-Voronoi tessellation and Pi -Lagrange inter- 
polation. The linear system obtained in such a way is 
solved by using the Cholesky method, from the Finite 
Element Library MODULEF [|o|. Examples of meshes 
with a 6*^ generation boundary are shown in Fig. ^. 
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FIG. 2. The building process of random von Koch curves. 
The same random process can create various interface to- 
pographies. They share the same size, the same perimeter, 
and the same fractal dimension. 




FIG. 3. A finite element mesh for the 6 generation von 
Koch random electrode. Top: an example of a finite element 
mesh for the 6*^ generation interface. Bottom: local zoom of 
the mesh. 

Computations were carried out for the two determinis- 
tic boundary geometries and the two random geometries 
of generation 6 shown in Fig. |^. The figure presents the 
isocurves of the potential for A = 0. Since the current 
density is proportional to the gradient of the potential, 
one can detect regions of large current density from the 
distance between two consecutive isocurves: the closer 
the equipotentials, the larger the current density. As ex- 
pected, most of the current flow through the irregular 
interface at the tips. This gives a very different current 
map for each geometry. Therefore, for the different elec- 
trodes the active zones are very different. 
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FIG. 4. Isopotential curves for von Koch deterministic and 
random electrodes with A = (Dirichlet boundary condition). 
The equipotential hues are the hues separating regions of ex- 
ponentially decreasing potential: V=l at the bottom then 
1/2, 1/4, 1/8,... The current density is proportional to the 
gradient of the potential. The current is then large in re- 
gions where the curves are close. Note that the current flows 
through the interface primarily at the tips. These active zones 
are found at very different locations for different electrodes. 

The second type of electrodes to be compared is shown 
m Fig. |. The top fi gure shows the second generation 
of a deterministic fractal electrode with dimension Df = 
In 16/ In 8 = 4/3 while the bottom represents a particular 
self-avoiding walk with the same 4/3 fractal dimension. 
Both electrodes have the same perimeter and the same 
smaller cut-off. Here, even more than above, the active 
zones are totally different. 

For each geometry, the impedances have been com- 
puted for an extended range of the surface resistivity r. 
The results are shown in Fig. ^ for two categories of ge- 
ometries : 6*^ generation of von Koch electrodes and 
the two electrodes of Fig. ^. The parameter A/ 1 = err /I 
ranges between 1 and 10^ for generation 6 and between 
10~^ and 5.10^ for the second type. The limitation of the 
range is due to limitations in computer time and memory. 





FIG. 5. Isopotential curves for deterministic and random 
electrodes of fractal dimension 4/3 with A = (Dirichlet 
boundary condition). Same color code as Fig. ^ The ac- 
tive zones are entirely different. 
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FIG. 6. Top: Plots of the electrode impedance Zs 
a function of A// = tct//, for various deterministic and ran- 
dom geometries. Note curve similarities despite very different 
current maps. Bottom: Plots of the "equivalent length" of 
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It is striking that, despite very different current distri- 
bution in the bulk and at the interface, the impedances 
are very close for all values of the surface resistivity. The 
behavior of different interfaces are nearly indistinguish- 
able: random and deterministic interfaces behave in the 
same manner. This could be considered as a partial an- 
swer to the question " Can One Hear the Shape of an Elec- 
trode?", addressed in |§Jl^. In this frame, the main pa- 
rameters drawn from practical impedance spectroscopy 
measurements would only be the size, the perimeter and 
the equivalent fractal dimension of the interface. 

A more demanding comparison between the 
impedances can be made by comparing the values of r/Z 
as shown in Fig. |[ This quantity can be identified as an 
equivalent active length Leg. [|l^. One finds three succes- 
sives regimes, A <l; l< A < Lp^ and finally Lp< A sepa- 
rated by smooth crossovers. These regimes can easily be 
compared to the so-called " land surveyor approximation" 
[|l3| . This method allows one to compute Zspec. through 
a finite size renormalization of the interface geometry, 
without solving the Laplace equation. For small r (or 
yl << i) there is a linear regime in which Zspec. is pro- 
portional to r, that is Zspec. = '^/Leq. with Leq. ~ L . 

For values of A> I there is a fractal regime in which, in 
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first approximation, Zs^ec. 

= {r/A){l/L){A/l)^/^f and 
Leg. = L{A/l)^^f~^^^^f (for more detailed expressions of 
tiie exponents, see |]l^-|l6[). Finally, for values of A much 
larger than the perimeter length L^, the exact value is 
Zspec. = ^/Lp and Lgg. = Lp. These three asymptotic 
behaviors are shown in the figure and are found to match 
the numerical results with good accuracy. 

Note that the electrodes of Fig. ^ are in some sense 
"poor" fractals because the range of geometrical scaling 
is relatively small and it has been a matter of debate 
recently whether the fractal concept should be of any 
use when the scaling range of the geometry is too small. 
For the phenomena considered here, one can observe that 
the fractal description of this limited range geometry is 
really useful. 

In summary, one has shown on several examples that 
the net transfer across an irregular surface is nearly in- 
dependent of the randomness of its geometry, although 
it depends strongly on the geometry through its fractal 
dimension. The fact that the overall response remains 
the same indicates that, buried in the fractal descrip- 
tion, there exist the geometrical correlations that govern 
the overall effect of screening at different scales. In that 
sense the response is "universal" within a very good ap- 
proximation for the category of curves considered here. 
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